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In this paper we are discussing the question how a continuous quantum system can 
be simulated by mean field fluctuations of a finite number of qubits. On the kinc- 
matical side this leads to a convergence result which states that appropriately chosen 
fluctuation operators converge in a certain weak sense (i.e. we are comparing expec- 
tation values) to canonical position and momentum Q, P of one-degree of freedom, 
continuous quantum system. This result is substantially stronger than existing meth- 
ods which rely either on central limit theorem arguments (and are therefore restricted 
to the Gaussian world) or are valid only if the states of the ensembles are close to the 
"fully polarized" state. Dynamically this relationship keeps perfectly intact (at least 
for small times) as long as the continuous system evolves according to a quadratic 
Hamiltonian. In other words we can approximate the corresponding (Heisenberg pic- 
ture) time evolution of the canonical operators Q, P up to arbitrary accuracy by the 
appropriately chosen time evolution of fluctuation operators of the finite systems. 

1. Introduction 

One of the most exciting, recent experiments in quantum information science is the realization 
of quantum memory. In [I] it is shown that the state of a laser field can be stored in collective 
degrees of freedom of an ensemble of about 10 23 earth-alcaline atoms at room temperature. When 
the light is released after a few microseconds it is in a state which has a very good fidelity with 
respect to the original state. 

A crucial part of the theoretical description of these experiments is the Holstein-Primakoff 
transformation [2], which basically implies that under certain conditions the ensemble can be 
treated like a continuous quantum system interacting harmonically with the light field. Storing 
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and releasing of the light is then - roughly speaking - just a swap of two coupled oscillators; cf 
[3J for details. 

One of the most important assumptions we have to make for this approximation to be valid is 
invariance of the interaction under permutations of the atoms (i.e. the intensity of the laser field is 
constant over the region where the ensemble is spread) , and this brings the whole model in contact 
with mean field theory and mean field fluctuations. The Holstein-Primakoff transformation can 
be regarded in this context as a special instance of the fact that fluctuation operators (which 
measure the fluctuations of a mean field variable around its expectation value) satisfy in the 
infinite particle limit and under quite general conditions canonical commutation relations [3] . 

This observation provides a big motivation to revisit the subject of mean field fluctuations and 
to look how it can be used to describe the type of matter-light interactions used in quantum 
memory experiments. Most of the results available in the literature arc mainly targeted towards 
the needs of statistical mechanics [5] and are strongly connected to the realm of the central limit 
theorem [5]. This implies in particular that in the infinite particle limit we only get Gaussian 
states. For the type of applications we have in mind this is a too severe restriction. To remove it 
is one main goal of this paper. 

To this end we are going to ask how a continuous quantum system can be simulated by mean 
field fluctuations of a finite number of qubits. This question involves kinematical as well as dy- 
namical aspects. On the kinematical side this means in particular to give the above statement 
about commutation relations of fluctuation operators a more detailed and mathematical precise 
meaning. This is done by a convergence result (Theorem 12. 9|) which states (roughly speaking) 
that appropriately chosen fluctuation operators converge in a certain weak sense (i.e. we are com- 
paring expectation values) to canonical position and momentum Q,P of one-degree of freedom, 
continuous quantum system. This result is substantially stronger than existing methods which 
rely either on central limit theorem arguments (and are therefore restricted to the Gaussian 
world) or arc valid only if the states of the ensembles are close to the "fully polarized" state. 
As a side result we will see that the mixedness of the one particle states (i.e. the restriction of 
the global state to a single particle) is related to an effective H in the corresponding continuous 
system. 

Dynamically this relationship keeps perfectly intact (at least for small times) as long as the 
continuous system evolves according to a quadratic Hamiltonian. In other words we can approx- 
imate the corresponding (Hcisenbcrg picture) time evolution of the canonical operators Q,P up 
to arbitrary accuracy by the appropriately chosen time evolution of fluctuation operators of the 
finite systems (Theorem l2.10[) . This is, again, a big improvement over existing results where only 
special Hamiltonians and weaker types of approximations are discussed. 

Before we start to present the details, let us have a short look on the organization of the 
paper. Section [2] collects all main results of the paper without proofs and keeping the technical 
overhead as small as possible. It consists of three subsections starting with [23] where fluctuation 
operators are reviewed and slightly reformulated for the needs of this paper. Subsection 12.21 
discusses Schwartz operators recently introduced in [7] . This class encompasses all states of the 
continuous system which can be simulated with our method. In Subsection 12 . 31 we present a new 
representation of permutation invariant states of many qubit systems and introduce a special type 
of convergence. It plays a crucial role in the two main theorem which are stated and discussed in 
Subsection 12.41 The rest of the paper is devoted to proofs (Section [3]) and an outlook to further 
research (Section!?]). 
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2. Basic definitions and main results 
2.1. Fluctuation operators 

Let us start to setup some notations and to recall a few known results which will be important 
for the rest of this paper. Therefore lets denote the space of (bounded) operators on the Hilbcrt 
space % — C 2 by B(H). For each A g B(H) and a fixed density operator 9 (the reference state) 
on H define the fluctuation operator Fm(A) £ B(H® M ) by 

1 M 

F M (A) = -= y £ / {A^-'Tr(Aff)l) (1) 



where 

A(*) = i»(*-i) 8 ^ 8 i^-*. (2) 

The Fm{A) are important in mean field theory and describe the quantum fluctuations around 
a mean field observables; cf [5] for details. For us the Fm(A) will serve as approximations (in a 
sense we will make clear below) of canonical operators of a continuous quantum system. To this 
end the central definition is the following: 

Definition 2.1. A sequence of density operators pm, M £ N has \/^-fluctuations if 

lim Tt(pmFm(A 1 )---F m (A k )) (3) 

exists and is finite for all K € N and all A\, . . . , Ak £ B{l-L). 

For the rest of this paper let us choose a reference state 9 satisfying (where er Q , a = 1,2,3 
denote Pauli matrices) 

Tr(0 CTl ) = Tr(&7 2 ) = 0, Tr(0a 3 ) = 2A, A e [0, 1] (4) 
and consider fluctuation operators corresponding to o\ and 02 • 

Qm = Fm(^]=^= PM = Fu(^)=±m (5) 



here Lm,ol denote global angular momentum (or global spin) operators given by 

W/ = jE' T i'' a = 1,2, 3. (6) 
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The operators Qm and Pm satisfy the commutation relations 



[Qm,Pm] - = -^=F N (a 3 ). (7) 
Z\f ivi 



Taking expectation values and the limit M — > 00 on both sides leads to 



lim Tr( pm[Qm, Pm}) - i\= — 7= lim Tr(p M F M (0-3)) = 0. (8) 



This indicates that Equation ([7]) converges in a certain formal way towards the canonical commu- 
tation relations. A more precise version of this statement is the following (the proof is postponed 
to the appendix): 
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Proposition 2.2. Consider a sequence pm £ B{H.® M ), M £ N of density operators with \fM 
fluctuations. Then there is a Hilbert space Jioa , & density operator poo and two symmetric opera- 
tors Qoo,Poa with the common, invariant, dense domain D such that for any polynomial f{q,p) 
in two non-commuting variables q,p the following statements hold: 

1. Poo D c D, 

2. The closures of f{Q 00 ,P 00 )p 00 and p 00 f{Q 0Q ,P 0Q ) are of trace class. 

3. ]im M ^ooTr{f(QM,PM)PM) = T^(f(Qoo,Poo)Poo) =r fr{Poof(Qao,Pao)), 
4- [Qoo,Poo]0 = iA0V(/)e£. 

Condition [1] and [5] arc of technical nature and needed to guarantee that [3] is a mathematically 
well defined expression. We come back to this point in the next subsection. Item[3]itself establishes 
the operators Qoo and as a form of weak limit of the fluctuators Qm,Pm, and condition 
U states that these limit operators satisfy the canonical commutation relations with A from 
Equation ^ as an effective h. We are tempted to interpret this result by saying that "the 
collective fluctuations of an ensemble of two-level atoms behave in the infinite particle limit like 
a continuous quantum system with one degree of freedom." Unfortunately this interpretation is 
premature, since condition [3] is not sufficient to conclude that Qoo and Poo are canonical position 
and momentum. It is not even guaranteed that Qoo and Poo are essentially self-adjoint on the 
domain D. Hence the crucial question is: 

Question 2.3. For which sequences pm £ B(H® M ), M £ N Provosition \2.2\ holds with 

1. "Hoo = L 2 (R) - i.e. the Hilbert space of square integrable functions on R, 

2. D = <S(R) - i.e. the space of Schwartz functions on R, 

3- Qoo = Q with (Q(j))(x) = x(f>(x) € <S(R) - i.e. Q is the canonical position operator with 
<S(R) as its domain. 

4- Poo = AP with P<fi = —i<p' V</> S <S(R) - i.e. P is the canonical momentum operator with 
<S(R) as its domain. 

One way to answer this question is to be more restrictive with the selection of the sequences 
of states pm, M £ N. An important case arises if take a translationally invariant state uj of an 
infinite spin-chain and define the pm in terms of the restrictions of u to a sub chain of length 
M. It is shown [S] that in the case where u> is exponentially clustering (i.e. correlations decay 
exponentially fast), Proposition 12.21 holds with Qoo and Poo being the canonical position and 
momentum operators and poo being a Gaussian state. 

In this paper, we want to discuss a slightly different approach which puts more emphasis on the 
non-Gaussian world. In the following two sections we will present in particular a different form 
of convergence which is tighter than pointwise convergence (as used in Definition 12. ip . and which 
can therefore guarantee that the M — > oo limit always leads to the Schrodinger representation 
of the CCR. 

Before we come to this point let us add a short remark about the role of the reference state 
6. To this end let us calculate first the expectation value of one fluctuation operator Fm(A). For 
an arbitrary pm we have 

Tr(^W) = Tr(^ M ,,) (9) 
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where 8m,i is the restriction of pM to the i th site (for a permutation invariant pu ah &M,i 
coincide). Hence 

1 M 

Tr(F M (A)p M ) = y/MTr((e M - 9)A), 9 M = jjY, 9m ^ ( 10 ) 

i=l 

Hence, liniAf->oc Tr (Fm (A) pm) < oo for all A £ BifH) can hold only if 

lim M = 6 (11) 

M — >co 

is satisfied. In other words: if a sequence pu has y/M fluctuations, its averaged one-site restriction 
in the limit M — > oo is given by 0. With the choice from Equation Q this implies that the 
parameter A - which plays the role of an effective H in Condition 2] of Proposition 12.21 - is a 
property of the sequence, while the operators Qmi Pm do not depend on A. 



2.2. Schwartz operators 

Conditions [1] and [2] of Proposition 12.21 are needed to guarantee that the traces on the right hand 
site of the expression in item [3] are well defined. They are, however, quite restrictive and lead to 
a new class of operators which are called Schwartz operators. 

Definition 2.4. An operator p € B(%oo) is called a Schwartz operator if P a Q l3 pQ 13 P a is for 
all a, a' ,j3,l3' £ No a trace-class operator. The set of all Schwartz operators will be denoted by 

SiUoo). 

Schwartz operators are introduced and discussed in [7J. For our purpose only few properties 
are needed which we will present in the following. For more details (in particular and the many 
degrees of freedom version) we will refer the reader to |7J . Most important for us is the fact that 
<S(%oo) admits a natural topology. 

Proposition 2.5. For each a, a' ,f3,(3' G No the functional 

SiHoo) 3 P ^\\p |Uw,H| P a Q' 3 pQ f3 'P a ' ||iG M+ (12) 
is a seminorm. 5(Hoo) together with this family is a Frechet space. 

In passing let us add the remark that we could replace the trace norm in (|12j) by any other 
p-norm (with p > 1 including the operator norm) and we would still get valid families of semi- 
norms (including in particular the property that all the scminorms are finite for all Schwartz 
operators) defining the same topology as the family we have chosen; cf [7] for details. 

There are a number of alternative characterizations of Schwartz operators and some of them 
we will encounter later in this paper. For now we only need a statement about some matrix 
elements. To this end let us introduce some notation first. The usual creation and annihilation 
operators are 

a=-j=(Q + iP), a t = _L(Q_ i p) (13) 

both with S(M.) as their domain. They give rise to the number operatoiQ 

N = ata, Nijj n =nip n , n e N Q (14) 

and its eigenfunction ip n , n € No, i.e. tp n denotes the n th order Hermite function. Now we have 
the following proposition [7J . 

1 The overline means - as usual — the closure. 
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Proposition 2.6. A trace class operator p is a Schwartz operator iff 

sup \(xp n , fnl> m )\ ■ (M + \m\) k < oo VfceNo. (15) 

This implies in particular that all operators p where only a finite number of matrix elements 
(ipm pipm) are different from 0, are Schwartz operators. 

2.3. Permutation invariant states 

Let us come back now to Question 12.31 The purpose of this section is to define a notion of con- 
vergence which guarantees that the operators Qoo and Poo in Proposition 12.21 can be chosen as 
canonical position and momentum Q and XP. To this end let us have a another look at Defi- 
nition [2TT] The crucial objects in Equation J2J are expectation values of products of fluctuation 
operators Fm(Am)- The Fm(A), however, commute with all permutation operators. Therefore 
the expectation values of the form Tr(i<M(^4i) ■ ■ ■ Fm(Ak)pm) do not change, if we replace the 
density operator pm £ £>»('H® M ) with its average over all permutations. In other words: As long 
as we are only interested in expectation values we can assume without loss of generality that pu 
is permutation invariant. 

The general structure of permutation invariant density matrices can be easily deduced from 
representation theory of SU(2) and the permutation group (this is very well known stuff; cf. [H] 
for a review and further references). The M-fold Hilbert space %® M can be decomposed such 
that for any U £ SU(2) we have 

H® M = Uj ® K M j, U® M = TTj (U) <g> 1 (16) 

j 

where Uj is a 2j + 1 dimensional Hilbert space, ttj : SU(2) — > B(Hj) denotes the irreducible 
spin-j representation, Kmj is a multiplicity space carrying an irreducible representation of the 
permutation group Sm, and the index j runs over the integers 0, 1, . . . M/2 (if M is even) or the 
half-integers 1/2, 3/2, . . . , M/2 (if M is odd). A permutation invariant density matrix has the 
form 

0"-"'(/'^ di^T") (17) 

where pu,j is a density matrix on Uj and the weights wmj are positive real numbers satisfying 



E 



WM,j = !■ (18) 



J 

J 



If ffiMj = for some j we set pu,j = 0, too. 

In a similar way we can decompose the angular momentum operators L a ^. If we writtj^l 

= TTj(a a /2) £ B(Hj) for angular momentum in the spin-j representation we get 

£a,iV = 04 J ' ) ®l- (19) 
3 

(i) (i) 

Of special importance as well are the eigenvectors ij)n of Lg 

L^^P = (j - n)^p , n = 0, . . . , 2j (20) 



2 In abuse of notation we use the symbol for the representation of the group and the corresponding Lie-algebra. 
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and the ladder operators L± = L± ± iL^ 

L ( i ] ^ = Vn(2j-n + l)^ ) _ 1 L®1>P = ~ n)(n + (21) 

Th next important points concerns the relation between the Hilbert spaces T-Lj and Hoo. The 
former is the representation space of the spin-j irreducible SU(2) representation and contains 
the distinguished basis ipn. The latter was introduced in Subsection 12 . 21 as L 2 (R) and we chose 
the Hermite functions tp n , n S N as a basis (cf. Equation (|14|) ). Now we embed all the T-Lj into 

(1) 

Hoo such that ip„ becomes for all n = 0, . . . , 2j the Hermite function ip n . In other words, we 
consider the embedding: 

2j 2j 
n=0 n=0 

In the following we will identify T-Lj with its image under this isometry. In this way all operators 
Pm become finite rank operators on the same infinite dimensional Hilbert space Woo- The same 
is true for the angular momentum operators L& \ a = 1, 2, 3 and lS± . 

Another crucial step is a rcscaling of the parameter j = 0, . . . , M/2 (or 1/2, . . . , M/2) by 

_2i [Mxl 
Xj ~ M ■' ' 2 ' 1 j 

Hence pu,j becomes PM.Xj and if we extend it continuously to the interval [0, 1] we get a contin- 
uous function 

R M : [0, 1]3I4 R M {X) G S(Hoo), Rm(Xj) = PM.y (24) 

Note that Rm is not uniquely defined by pm, however we can choose such a function for any 
Pm (e-g- by linear interpolation). Also note that the definition of these functions requires the 
embedding of the Hj into Hoo as described above. 

The advantage of this formulation just introduced is that sums over j can now be reformulated 
as integrals with respect to the probability measure pu given by 

1 

f(x)p, M (dx) = y^w M ,jf(xj). (25) 

j 

To see this consider the expectation value of a permutation invariant observable (like angular 
momentum introduced above) Am = Am, j <8> I in the state pM- Using the decomposition in 
Equation (TT7]) we get 

Ty(pmAm) = y~] wmj Tr(p M ,jA M ,j) = / r Tr(R M (x)A M {x))pM(dx) (26) 
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where we have used again a continuous function Am '■ [0, 1] — > B(Hoo) satisfying Am{xj) = Amj- 
This motivates the following definition: 

Definition 2.7. Given pm, permutation invariant. We call a continuous function Rm '■ [0,1] — > 
S(Hoo) an integral representation of pm, if Rm{2] /M) = pM,j> where pmj € B^ifhij) are the 
density operators in the decomposition from Equation &17\) . 

Note again that an integral representation always exist, but is never unique. The latter is not a 
problem, though, since integrals as in Equation (|26[) only depend on pm and not on the integral 
representation chosen. 
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Finally let us have a look at the convergence of sequences of permutation invariant density 
operators P m, M G N. We can split this question up into two pieces: convergence of the sequence 
of measures p Pt M and convergence of the sequence of functions Rm- A particular behavior is 
described by the next definition, which captures the situation that in the limit M — > oo only one 
value x = A G [0, 1] is relevant, while all other x does not matter. Or to state it in a different 
way: For large M the state pu describes a very sharp angular momentum L? with j centered 
around j\ = AM/2. For these j the corresponding pM,j is close to p^ (in the topology of S(Hoo) 
which measures in turn the moments). All other pj are unimportant. 

Definition 2.8. We say that a sequence pm G B*('W® M ) of permutation invariant density op- 
erators converges at A G [0,1] towards a state p^ € S{'H ao ) if each pm admits an integral 
representation Rm such that the following conditions hold: 

1. The sequence of probability measures pm converges weakly to the point measure to the point 
measure at A that is 

lim / f{x)p M {dx)= f{\) (27) 

A/->oo J 

for all continuous functions f : [0, 1] — > R 

2. The set 

{Rm(x)\M €N,xe [0,1]} c S(Hoo) (28) 

is bounded. 

3. There is a neighbourhood of I a [0, 1] of A and a continuous function R^ : I — > 5(Hoo) 
such that R 00 (X) = p^ and 

lim Rm(x) ~ Roo{x) uniformly on I (29) 

M— >oo 

holds in the topology of Sfflca). 

2.4. The main result 

Now we can come back to Question l2.3l The next theorem simply says that the convergence given 
in Definition 12.81 provides (at least) a sufficient condition (all proofs are postponed to Section [3]). 

Theorem 2.9. For any sequence pu g B*{T-L® M ), MgN of permutation invariant states con- 
verging towards a poo G S{Hoc) at X € (0, 1], and all polynomials f{q,p) in two non-commuting 
variables we have 

lim Tv{p M f(Q M , P M )) = Tr( Poo fiQ^, Poo)). (30) 

M— >oo 

with Qoo = ^f~\Q and = VAP. 

Note that convergence at A in the sense of Definition 12.81 does not imply \J~M fluctuations. 
The reason is that we are not controlling convergence of the moments of the measures pm- 
Therefore expectation values of products involving Fm(o~3) can still diverge. It is easy to fix this 
by considering a stronger type of convergence for the pu- However, this is not needed for getting 
the proper limit of the Qm, Pm- 

One way to interpret Theorem 12.91 is in terms of quantum simulation: If we are interested in 
the expected value of an observable X = f{Qoo, Poo) of a continuous quantum system in a state 
Poo, we can measure instead the observable Xm = /(Qa/, Pm) of an ensemble of M qubits in the 
state pm and the error we make that way can be arbitrarily small, provided M is big enough. 
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It is possible to simulate each density operator p,^ G Sffloo) that way. To see this consider for 
each half-integer j the projection 

Ej : Hoc -> span{V)„ | < n < 2j} (31) 

and a sequence jm G R satisfying j'm G {0, . . . , M/2} if M is even, jm G {1/2, ... , M/2} if M is 
odd, and liniM^oo ^Jm/M = A. Now we can define 

= E juPoo E jM ® 1 G B(W jV ® /C MJm ) C £(?^ M ) (32) 

and 

Pm = Tr ^^ if Tt(pm) ^ 0, (33) 

where we have used the identification of SU(2) representation space Hj M with the subspacc 
Ejm'Hoc introduced in Subsection 12.31 If Tr(pju) = holds (this can only happen for finitely 
many M) we can just choose an arbitrary density operator. It is obvious that the sequence just 
constructed converges towards Poo at A and therefore it simulates p^ in the sense described in 
the last paragraph. 

The parameter A can be freely chosen in this construction (by choosing the sequence jm) 
without affecting the simulability of poo . This seems to indicate that A does not play an important 
role, but this impression is wrong. At the end of Subsection 12.11 we have seen that A describes 
(asymptotically, i.e. in the limit M — > oo) the noise in the one-site restrictions of the pm, ranging 
from "fully polarized" for A = 1 and "fully depolarized" for A = 0. 

At the same time it plays in the commutation relations in Proposition ^. 2l the role of an effective 
h, and this also holds for the operators Qoo and P^ in Theorem 12.91 However, in contrast to 
our expectation in Question 12.31 Qoo and are not exactly canonical position and momentum 
(except in the case A = 1), they differ by a factor A" 1 / 2 and A 1 / 2 respectively. Hence, although 
A has some properties of h the limit A — > oo is not the classical limit (unless we rescale the finite 
dimensional observables Qm and Pm appropriately). Instead all observables vanish if A — > oo. 
Nevertheless, the case A = is interesting (in particular because it is important experimentally 
[0]) and therefore deserves a more detailed study in a forthcoming paper. We come back to this 
point in Section 3J 

Now let us come back to the simulation point of view. We have argued that we can simu- 
late states, but what about time-evolutions? To answer this question, consider a self-adjoinl[^| 
polynomial h(q,p) and define the Hamiltonians Hm = /i(Qm,-Pm) and H = h(Q,P). Theorem 
12.91 shows that we can simulate expectation values of H by expectation values of Hm- But what 
about time-evolutions? This question is much more difficult. The first big obstacle arises from the 
fact that H (with <S(R) as its domain) is only symmetric but in general not essentially selfadjoint. 
It is even possible that no self-adjoint extensions exist at all [TU]. In these cases the unitaries 

U M .t = cxp(-itH M ) (34) 

(which always exist, since the underlying Hilbert space is finite dimensional) do approximate a 
reasonable time-evolution of the continuous system. However, what about the case where H is 
essentially self-adjoint on <S(R)? In this case we can define the weakly continuous one-parameter 
group 

U x .t = cxp(-itXH) (35) 

3 For each polynomial f(q,p) in two non-commutative variable we can construct the adjoint we apply complex 
conjugation to the coefficients and reverse the ordering of each monomial. A polynomial is self-adjoint, if it 
coincides with its adjoint. 
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and ask whether the finite dimensional dynamics Umj approximates the U\ t t with a certain 
error for a finite amount of time. The next theorem shows that this indeed the case for quadratic 
Hamiltonians. 

Theorem 2.10. Consider: A sequence pm £ #*(%® M ), MeN of permutation invariant states 
converging towards a £ 5("Hoo) o,t A € (0,1] , a second order, self adjoint, homogeneous 
polynomial h(q,p) and the corresponding unitaries, as defined in Equations fl5^| ) and h35}). Then 
there is to > such that for all t £ R with \t\ < to o,nd all polynomials f(q,p) we have 

lim Tr(U M ,tPMU* M J(Q M , P M )) = Tr(U x , tPoo U* x ?„)). (36) 

M— >oo 

im't/i Qoo = yXQ and P^ = s/\P. 

This theorem implies that the time evolution of the expectation values (f(Q, P)) of any ob- 
servable f(Q,P) can be approximated for short times and with a bounded error e by the time 
evolution of the corresponding expectation value {f{Qu, Pm)) of the finite dimensional system. 
The size of the error can be as small as possible provided the number M of qubits is big enough. 
This strengthens the simulation point of view introduced above (at least as long as quadratic 
Hamiltonians are considered). 

There is, however, one small problem with this result, and this is the time limit to- As long as 
we are considering an ensemble of fixed size M it is natural that we can keep a bounded error 
e only for a finite amount of time. The bound to in theorem, however, holds also in the infinite 
particle limit. This is counter intuitive since we would assume that we can always improve the 
quality the simulation by increasing the number M of particles. Therefore, it is likely that the 
restriction of a finite to is only a restriction of the proof and not of the result. 



3. Fluctuator dynamics 

The purpose of this section is to prove Theorem 12.101 and in this context we will learn more 
details about the dynamics of fluctuation operators. As another side effect, we will get a proof 
of Theorem l2.9l as well, since it arises from l2.10"l in the special case t = 0. Let us start with some 
additional notations. 

• We need the functions 

P : N 2 x [0, 1] -> [0, 1] : (M, n, x) ^ u (x, n) = 9 x (x - ^) 

with 6i{x) = ^xx[ ,i](x), (37) 

where xs is the characteristic function of the set S. Some of the arguments of j3 will be 
often omitted when confusion can be avoided. 

• Furthermore we introduce am '■ [0, 1] — > <S(%oo)>VM £ N by 

a M (x) = f3 M (x,N)a = af3 M (x,-N-l) (38) 
a* M (x) = /3 M ,A(z,N-l)a* =a*/3 M (x,N) (39) 

where a, a* are the standard annihilation and creation operators, respectively, and N = a* a 
is the number operator (cf. Scction [2.2p and f3 with an operator argument is understood in 
the sense of functional calculus. 
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We will also need the integrated version of 0^(1), a* M (x), which are given by 

aM = ~7TT = ~m^ M + iPM ^ a *M = ~7Tt = 4^( < 5 M ~ iF>M ^ ( 40 ) 

with 

1 M 

L M ,± = L M ,i ± iL M ,2 = g ^ a ± ( 41 ) 

j'=i 

with cr± = CTi ± ifT2 in terms of Pauli matrices, bjk, are related to the functions am(x), 
a* M (x) by 

a M = ® 1, = a^z) <g> 1. (42) 

i j 
where ajj is given as in Equation (|23[) by = 2j/M. 

The quadratic Hamiltonian iJ will be written as 

3 

H = ^ c kA k with A = a 2 , A x = aa*,A 2 = a*a, A 3 = a* 2 (43) 

fe=0 

and the finite dimensional version Hm becomes similarly 
3 

Hm = CkA M>k with A M ,o = a 2 M , A M ,i = a M a* M , A M .i = a* M a Ml A M<3 = a*^. (44) 

fe=0 

In analogy to Equation (|42p we can rewrite Hm as a direct sum 

H M = ($H M {x j )®l (45) 
3 

where 

3 3 

Hm (a:) = £ «\«w,, (a:)^ = £ c 3 A M ., (x), (46) 

with 

B M ,o(x) = Pm(x, N)f3 M (x, N + 1) B M<1 {x) = M {x, N) 2 

(47) 

B M , 2 {x) =/3 m (x,N-1) 2 B Mt3 (x) = p M (x,-N - l)p M (x,N -21) 

Finally we can write the unitary J/m,* from Equation (134j) in exactly the same way 

^M,f = 0Efo,t(aO®l (48) 

with 

t/ M ,iO) = exp(-iiff M (a;)). (49) 

Furthermore we define the multiindcx R £ { — 1, which enables the compact notation 

„ii — /1 (Ri)_(^2) (R d ) „(i)_„* „(-!)__ ^n^i 
a = a M a M ... a M a M - a M a M - a M (5U) 

and similarly for the quantities a £ S(Hoo). In addition we define \R\ = d and w(R) = 

Et Ri- 



ii 



• To enable compact notations in some proofs we will write objects describing the continuous 
system sometimes with the subscript oo, i.e. 

aoo(x) = y/xa, a^(x) = \/xa* , H (X {x)=xH, Uoo^ix) = U Xjt , (51) 
where x is as usual in the interval [0,1]. 
To prove Theorem 12. 101 we have to show that for each e > we can find an M e such that 

|TV(C/ M ,tPM^ jt O - Tr(U x , tPoc Ul t a R )\ < e (52) 

holds for all M > M e . Note that we have replaced the polynomials /(Qm, Pm) and f(Q, P) from 
Theorem 12. 101 by monomials a R j and a R with an arbitrary multiindex R. The general statement 
easily follows by linearity. To get such an estimate note first that we can rewrite traces of the 
form ^{Um^PmUm^) as 

Tr {UmjPmUmjCIm) = ^ w M (j) r TT(U M ,t( x 3) R M(xj)Ul Itt (x j )a M (x j )) (53) 

3 

= / Ti{UM.t{ x )RM{x)U* Mt (x)a M {x))pM{dx) (54) 
Jo 

where Rm is an integral representation of pu and \xm is the measure define in Equation (|25[) . 
With this we can rewrite the right hand side of ([5^]) as 

|Tr(C/ M>t PME/V4) -Tr(C/ x , tPoo [/ A >«(A))| < 

/ \TY{UM.t{x)PM{x)Ul lt {x){a B u (x) - a R (x)))\ f i M (dx)+ 
Jo 

/ \Tr({Ux,tPooUl t - t/ A /, f (a;)p 00 t r ; / . t (x))a^(x))|/i M (dx) (55) 
Jo 

Note that we have used here a^o(A) and aoo(A) from Equation (|51|) . They are creation and 
annihilation operator belonging to Qoa and P^ from Theorem 12.91 and 12.101 

The discussion of the two integrals on the right hand side is now broken up into four steps. 

• In Subsection 13.11 we will collect additional material about Schwartz operators, which is 
used throughout this section. 

• Then we show that the sequences t a R j{x)l) 'm ,t^n converge for any Hermite function to 
U* t a R U x ^ n ] cf. Subsection^ 

• In Section [3751 we look at operators (N + V) v ^^m t a Ai( x )UM,t and show that pgN can be 
chosen such that the sequence of their norms is uniformly bounded in x. This allows us to 
trace convergence of unbounded operators back to convergence of bounded operators; cf. 
Subsection 13.31 

• The two previous steps together allows us to bound the integrands of the integrals in 
Equation (|55|) and this leads to the required estimate; cf. Subsection 13.41 



12 



3.1. More on Schwartz operators 

Let us collect some properties of Schwartz operators which are needed for the proof. The standard 
reference is [JJ. The first statement says that the space S(1-Loo) is stable under multiplication with 
polynomials in Q and P. 

Proposition 3.1. Consider a polynomial f(Q,P) in position Q and momentum P and a 
Schwartz operator p. The products f(Q, P)p and pf(Q,P) are closeable operators and their clo- 



sures are again Schwartz operator^]. The two maps 

SiHoo) 3 P ^ f(Q,P) P eS('H 00 ) (56) 
S{Hoo)3 p^ P f{Q,P)eS(H 00 ) (57) 

defined that way are continuous. 

Proof. This is an easy consequence of the definition. □ 



Furthermore the space S{'H 00 ) is closed under multiplication. In other words [7] 

Proposition 3.2. The product of two Schwartz operators is again a Schwartz operator. 

The next result is an alternative characterization of Schwartz operators. To this end let us first 
introduce for any x = (x±, X2) € K 2 the Weyl operator acting on tp € L 2 (]R) by 

(W(xi,a?aM(x) = e- ix ^ 2 / 2 e iX2X tp(x - x x ). (58) 

The W are unitary and for any density operators p they give rise to the characteristic function 

R 2 9i4 p{x) = Tr(pW(x)) G C (59) 

Now we have [7J: 

Proposition 3.3. A density operator p is a Schwartz operator iff its characteristic function p 
is a Schwartz function. 

Now let us come back to the quadratic Hamiltonians H and Hm{x) defined in Equation ([^3"]) 
and (j^r?)) . and the corresponding time evolutions U\ t t and Uu,t(x) from (|35[) and (|%9"]L 

Lemma 3.4. For each Schwartz operator p G S(TLoo) and each multiindex R G { — the 
following statements hold: 

1. U\,tpU^ t and UM,t(x)pUl I t (x) are Schwartz operators. 

2. t a R U\.tp and t (x)a R UM ,t(x) p are Schwartz operators (and in particular trace class). 

3. We have 

Tr(Ul t a R U x , t p) = Tr (a H U x , tP Ul t ) (60) 
Tv(U* M /x)a R U A u(x)p) =Tr(a R U M ^x)pU* M j(x)). (61) 



4 In slight abuse of language we will shorten in the future statements of this form as: "... is a Schwartz operator". 
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Proof. To prove the first statement in [T] note first that for quadratic Hamiltonians the canonical 
operators Q t = J7 t * x QU t ,\ and P t = £/ t * x PU t ,\ evolve according to the corresponding classical 
equation of motion (this is well known or otherwise an easy exercise). This implies for the Weyl 
operator that U% x W(x)Ut,\ = W(F t x) holds with a linear map F t (the classical phase space 
flow). Hence we get 

(uZpU* t ) = poF t , (62) 

which implies that (U X jpU x t )is a Schwartz function and therefore U Xj tpU x t is a Schwartz oper- 
ator by Proposition 13.31 

Now consider the unitary UM,t(x). Since its generator Hm{x) satisfies Ej x HM{x)Ej x = Hm(x) 
with j x = M^il f r om Equation (|23[) and the projection Ej x defined in (|3"Tj) . we can rewrite U M,t{x) 
as 

U M ,t{x) = (1 - E Jx ) + U M , t (x), U M ,t( x ) = E i*U M ,t( x ) E 3*- ( 6 3) 

Hence 

U M ,t(x)pUl It (x) = 

(1 - E jx )p(l - E jm ) + (1 - E Jic )pill ut {x) + UmAx) p (1 - E 3 J + U M .tpU* M ^ (64) 

Since p G iS('Hoo) by assumption, Proposition 12 . 61 immediately implies that (1— £ , Ja .)p(l — Ej x ), 
(1 — Ej x )p and p(l — Ej x ) are all Schwartz operators. The same is true for Umj an d t ( a g a i n 
by Proposition 12. 6p . Hence by Proposition 13.21 all operators on the right hand side of (|64[) are 
Schwartz operators and therefore UM,t{x)pUlf t {x) is a Schwartz operator as well. 

To show statement [3] note first that a R Ut t \pU^ x is a Schwartz operator - this follows from 
Proposition 13. II and statement [T] above. Applying item [1] again we see that 

Ul t a R U x , tP = Ul t (a R U t , xP U t * x )U t , x (65) 

is a Schwartz operator as well. Similar reasoning holds for t {x)a R l) M,t(x)p. 

The last statement is an easy consequence of item[5]and elementary properties of the trace. □ 

3.2. Hermite functions and analytic vectors for quadratic Hamiltonians 

The subject of this section is convergence of sequences a UM,t(x)ip n , M € N to a s U x .tip n - This 
will lead to a first partial result in Proposition l3. 101 stating that the UM.t(x) converge strongly to 
U x .t- The main technical tool will be the fact that quadratic Hamiltonians admit (finite) linear 
combinations of Hermite functions as analytic vectors (cf. [TTJ Ch. 5]). The first step is a general 
estimate. 

Lemma 3.5. Consider a Hermite {unction ip n and a multiindex S G { — l,l} 2d . Then the fol- 
lowing bound holds for all M G N U {oo} and x G [0, 1]. 

\\a s 'H M (x)ip n \\ < 2 3d+ td!m!(32max| Cj -|) m (66) 

with the maximum max \cj \ taken over the j = 0, . . . , 3. 

Proof. It is sufficient to show Equation (|66|) for the Hamiltonian H = H^l) instead of Hm(x) 
since all bounds hold also when we substitute aM(x) R for a R . The former differs from the latter 
only by a prefactor between and 1. Now consider the bound 

[\\a R tp„\\ < 2%4 d d\ . (67) 
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\\a R ip n \\ < y/(n + l)(n + 2) ■ ■ ■ (n + 2d) < V '2 n + 2d ^/(2d!) < 2t+ d 2 d d! = 2^4 d d! , (68) 
where to see the first inequality we recall 



a*ip n = y/n+lipn+i, aip n = vWVs-i , (69) 

from which the inequality is clear as it is an equality exactly in the worst case of a R = (a*) 2d . 
For the next two inequalities we used 

(P + 9)l < 2P+q (70 ) 
p\q\ 

(which holds since the lhs. is a term in the binomial expansion of the rhs.) with p = n, q = 2d 
and p = q = d, respectively. 

Let us now introduce the multiindcx j S {0, 1, 2, 3} m to enable the compact notation 

H m = CjAj with x 3 = x n x n . . . x jm (71) 
j 

Now, taking into account that 

| Cj | < | max Cj \ m Vj; (72) 

J=0,...,3 

and that Aj = a R for each j with an R satisfying \R\ = 2m, thus 

a s Aj = a R ' with \R'\ = 2{d + m) (73) 

we can calculate 

\\a s H M {x) m ^ n \\ < J2 |cj|||a 5 Aj^ n || (74) 
j 

< 4 m (max \cA) m 2%'i d+m (d + m)\ (75) 

< (32max(| Cj |)) m 2 3d+ *m!d! , (76) 

where in the second inequality we inserted the bound (|67|) for \R'\ = 2(m+d). The third inequality 
is due to the application of (fTO)) with p = m,q = d. □ 

Lemma 3.6. The bound in id 6} ) also holds if |s| = 2d — 1 

Proof. We obviously have 

\\aa s H M {x) m i> n \\ > \\a s H m 4> n \\ , (77) 
so the statement follows from Lemma l3~51 □ 
Lemma 3.7. For each M,n G N and x € [0, 1] the following bound holds 

\V^-^M(x,n)\<J^ (78) 
Proof. Recall the definition of j3 (cf. Equation ([3T)l : 

13 : N 2 x [0,1] [0,1] : (M,n,x) ^ (3 M (x,n) = 6 x (x - with 6 x {x) = y/x X[ o,i](x). (79) 
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If x < n/M the argument of the root is negative and we get /3m(x, n) = 0, hence 

\y/x- p M {x,n)\ = \fx< 
as stated. If x G [n/M, 1] we can check that the function 



I n/M, l]3xi-ty/x — \/x — — G 



(80) 



(81) 



is positive and monotonically decreasing (the latter can be easily seen by looking at the deriva- 
tive). Hence the biggest value is achieved at the beginning of the interval (i.e. x = n/M) and we 



get again the statement of the lemma. 

Lemma 3.8. For any pair of multiindices S, R and n G N we have 



M 



lim ||a s (<4(x) - x lRl/2 a R )iP n \\ = 



□ 



(82) 



uniformly in x. 

Proof. Note first that for r G {-1, 1} 



(3 M [x,n + 



r - 1 



r + 1 



-1pn+r 



The general case can be schematically written as 

d 

a s {a R i (x) - x W,2 a R )i> n ~ J[ y/n+ P] (f3 M (x, n + pj) - a S 7/W g 



(83) 



(84) 



with <i = |i?| and pi,p\,q G N can be easily written in terms of Wj(R) = y^ —1 r j with i? = 
(r"d,rd-i, . . . , ri). Note, that the omitted factor of proportionality is M-independent. Now we 
have according to Lemma 13.71 for any p G N 



77. + p 
M 



Hence lim^/^oo /3m{x, n + p) = *Jx uniformly, which proves the statement. 

Now, we can prove convergence of UM,t(x)tp n for small times. 
Lemma 3.9. Let x G [0, 1], n G N and \t\ < (32 max \c 3 . Then 



(85) 
□ 



A I 



lim \\a s Ut,Ai(x)->p n - a s U t , x ^n\\ = 



(86) 



uniformly in x. 



Proof. Note first that the expressions in Equation (|8"6")) are well defined, since V r M,t(x)ip n and 
Ux.t'i'n are Schwartz functions (this follows from Lemma l3.4j) and therefore in the domain of a . 
Now proceed with q = 32\t \ max \ cj\ < 1 which implies 



|a s ^PVn||<2 3d -td!g" 



(87) 
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The sum over mSNis then convergent since the rhs. of (|87|) is also (since q < 1). Hence we can 
write U x jipn in terms of the exponential series 



U x , t ipn = exp(-itxH)4' n = 2J 



(~itxH) m ^ r , 



Now note that U x ^ n € <S(R) and therefore a U Xjt if) n is well defined. Hence we get 



'Uxjipn = a s exp(-itxH)4> n = ^ 



a s (-itxH) m iP n 



(89) 



n=0 



The last equality holds (although a s is unbounded) because a s is closable and convergence of 
the series is known from (|57|) . 
Now, consider 

\\a s (U M ,t{x) ~ U x ,t)ip n \\ = || exp(itH M (x)) - exp(itxH)]ip n \\ < 



E 



a s (itH M (x)Y 



-</Y 



v-^ a s (itxH) m , 
i ^ 



E« 



s (itH M (x)) m — (itxH) 1 



■ipT, 



(90) 



Now, for all e > there is /i £ G N (independent of M and x) such that both the first two terms 
are smaller than e/3, since they are bounded by the remainder of the geometric series c.f. (|57|) 
and the remark after it. For the third term we write 



\^ s (itH M (x)) m -{itxH) m , 

2^ a m \ V* 



m</j, e 



< \\a S (H M (x) m -x m H m )^ 



\t\ 



^E E ||« 5 -s 1 * 1 ' 2 **)^ 

m<M E fi£{-l,l} 2m 



(91) 



(92) 



The last expression contains finite number of terms, each of which converges (uniformly in x) 
to zero when M — > oo due to Lemma (|3.8[) . hence for M large enough the last term in (j9"0|) is 
smaller than e/3 as well. □ 

This lemma leads to the following strong convergence result, which is a first step towards 
Theorem [2~TU1 



Proposition 3.10. For any fixed r\ £ T-L c 



lim \\(UmAx)-U x M\=0 

M-»oo 



(93) 



holds uniformly in x. 

Proof. Consider rj € Woo with ||?^|| = 1 and its expansion in terms of ip n : r\ = ^ n r)( n 'ijj n . For 
any e > we can decompose r\ = rjx + 772 such that 



E »7 (n tyn, ^ = E H^ll < e 



n<n e n>n € 



(94) 
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Now wc write 

\\(U M , t (x) - U Xit )r)\\ < \\(U M , t (x) - U X M\ + \\{U M ,t(x) - U x , t ) m \\ (95) 
< V {n) \\{U M ,t{x) - U m ,M\ + 2e (96) 

n<n e 

According to Lemma (|3.9[) for each n < n e there is an M n ^ such that for M > M„ jfr and all x 
we have \\(U M ,t( x ) ~ U Xtt )ip n \\ < e/n t Since ||r/|| = 1 we get 

\\(U M ,t{x) - U x . t ) m \\ < n e — = e VM > M t = max M„, 6 . (97) 

Together with (|96|) this leads to limM^oo UM,t(x)r/ = U x ,tf\- Since 77 g "Hoo was arbitrary, the 
statement follows. □ 

Now we are aiming at convergence (in a sense we will make precise later) of U m ,t{x) pU m ,t{x)* 
to U x .tpU* t . The next lemma is the first step. 

Lemma 3.11. Vn S N we have 

Hm U M ^x)a s U* Mt {x)^ n = U x>t a s Ul t i )n (98) 

uniformly in x. 

Proof. As in Lemma [3.91 note first that the expression in (f$g|) is well defined since t (x)ip n 

and Uxjtpn are m the domain of a s . 

From the same lemma we have in addition that $?M,n = aS U1^ t (x)ip n converges uniformly to 
$„ = a s U* t ip n . Hence 

\\U M>t (x)* M , n - U x , t <S> n \\ < \\U M ,t(®M,n - *«)|| + \\(U M .t - U Xtt )* n \\ < 2e (99) 

if M > max{Af e ,M e '} and for arbitrary x, where M e comes from the strong convergence of the 
second term proved in Prop. ROTH whereas M' e from the above mentioned convergence. □ 

3.3. The key estimate 

The purpose of this Section is to prove the following Lemma which will allow us to trace the con- 
vergence of sequences of unbounded operators t (x)a UM,t{x) back to convergence of bounded 
operators. 

Lemma 3.12. 3p,p' 6 NU {00} and K G M + such that there is a bounds 

\\(N+pl)-iU M ,t(x)a S UM !t (x)\\ < K (100) 

and 

\\U M Ax)a s U* M , t (x)(N +p'l)-^\\ < K' (101) 
for sufficiently small t and VM € N and for all x G [0, 1] . 
Proof. Since ip n , n £ N is a basis, it is sufficient to show that 

\\(N+pl)-%U M , t {x)a s U* AU (x)ip n \\ < K (102) 
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and 

\\U M ,t(x)a S U* ALt {x)(N +pI)-*Vn|| < K (103) 

hold. Furthermore, it is sufficient to show this for n > no since the overall bound will then be 
given by max{K, K„ \n < hq} where we denoted K the bound for n > no, and K n is a bound for 
fixed n's in (|102p . which exists due to Lemma r3.Hl For the proof of p02|) and (|103p Vn > no with 
a suitable no € N we proceed by introducing the spectral projections of the number operator 
E( a ,b) = X[a,b] (N) and proving some technical lemmas. 

Lemma 3.13. There is an n e £ N and t t G R + suc/i that 

||(l-£ [f) at])$M,n|| < e wit/i $ M ,n = a s Uli f1 .{x)ip n (104) 
VM e N, |t| < t e , n > n e . 

Proof. Since a H^(x) is a polynomial in aji/(i), a^(x) of order not greater than d + 2m (with 
d = \S\), we have 

(1 = if ^>d + 2u, (105) 



2 



where n denote the partial sums 



< - e -"cyv doe) 

m<jU 

Hence for the limit of the series the following formula holds: 

(1 - %, ¥] )<&M,„ = ((1 - E { , m )^ M . n - $it M)/4) ) • ( 107 ) 

Let us estimate 

ll(l-% >¥ ])(*^n-*Kr M)/4) )|| < ll(*M.„-4t" 2d)/4) )H ( 108 ) 

< £ ^MM, r (109) 



< 2 3d+ *d! ^ q m (110) 



i-2d 



with q = 32|t| max|cj|. With sufficiently small t, c.f. Lemma T3.91 that is, q < 1 is assumed, we 
can sum the geometric series explicitly: 

„-2d +1 3d o — d 1 

2 3d +f d!^ = 1 exp( logg) exp[n(log y/2+- logo)] = „ e^ 2 -*" (111) 

1 — g 1 — q 2 4 

where the constants Ki^ q can be read off from the formula. Now shortening the time interval 
such that q < 1/4 we find make K2,q < 0. Thus, we can conclude that there is an n^ q such that 
the statement of the Lemma is fulfilled. □ 

Lemma 3.14. 

d 

/ 3n \ 2 

\\E [§ ^ M>n \\<( — +2d) (112) 
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Proof. Let us expand 5?M,n an d UM,n(x)ipn in terms of Hermitc functions 

®M,n = C ^k, $ = ^2 Ck^k ■ (113) 
k k 

Then C k+w(S ) = Hk,S)c k holds with A(k,S) satisfying (k - d) d / 2 < A(k,S) < (k + d) d / 2 , 
where we recall | | = d. Since the relation — d < w(S) < d also holds, the condition k + w(S) € 
[n/2, 3n/2] implies fee [n/2 — d, 3n/2 + d}. From these we get 

||%. ¥] <f>M,„!l 2 <El c ^)l 2 (Y + 2d ) d ( 114 ) 

i 

and since J^ fc |c/c| 2 = 1 due to \\ip n \\ = 1 an d the unitarity of UM,t(x) we finally get to the 
statement of the Lemma. □ 

Now let us return to the bound in Equation (|102j) and decompose the norm in the following 
way. 

\\(N+pt)-SU M ,t*M, n \\ < W^+P^E^^Um^^mA 

+ \\(N + pl)-Hl-E [hS £ ] )U M , t E [ , m <Z> M , n \\ 

+ ||(N+pl)-ic/M, f (l-%,^])*A/.„|| (H5) 

The third term can be estimated as 

\\(N + pl)-%U M , t (l-E [hm )$ M>n \\ < ||(N + pl)-«||||(l-E [fi¥] )$ M ,„|| <p-h (116) 
whenever n > n e due to Lemma 13.131 For the second term we write 

+pl)-i{l- E^^UmjE^^mA < 

\\(N+pt)-%\\ max J|(l-^*.*l)^,*^llll B f*,*l*^,»ll < 



_ 3ti 

L 2 ' 2 



p- P -K^ q er K ^ + , (H7) 

where we used the Lemma f3. 141 and the technique of the proof of Lemma f3. 131 for the middle and 
right norms in the second line to arrive at the third. Note, that for k £ [n/2, 3n/2] we have that 
[fc/2,3fc/2] D [n/4,9n/4], which justifies using the same constants Ki q as in Lemma T3.13I It is 
now clear, that increasing n makes this bound arbitrary small. Finally, note that (N + pl)~ p / 2 
commutes with -E[ n /4.9„/4] , and since the latter is a projection we can write for the first term 

UN + piyiE^^UMjE^MfiuA < 

\\E [im {N+pt)-iE [im \\\\U M M E li,^MA > ( 118 ) 

Now, for the norms we have 

d 

||% )¥] (N+pl)-ts [fi¥] ||<2f(n + 4p)-«, ||B [f M]$M,n|| < (y + 2d ) 2 , ( 119 ) 

the multiplication of which, after substitution p = d, gives % d r d l 2 < 6 d since r = (n + 4e?/3) /(n + 
Ad) < 1. Hence we have found K = 6 d + 2e, which leads to ^TU7\i and therefore ([MI)) . To 
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get (|lU3p and (|1U1|) we can proceed in exactly the same way. The only difference is that the 
operator (N + p'l)~^~ is acting directly on ip n and not at UM,tE[2t,s^]^M,n- Therefore we only 
have to replace the first norm in (|119p by (n + p') p I 2 to get the desired result and the proof is 
complete. □ 

Now we are ready to prove a convergence statement for the time evolution UM,t{x)pU^ t (x) 
which brings a step closer to our goal. 

Proposition 3.15. For any Schwartz operator p G iS('Hoo) we have 

lim Tv[a s U M ,t{x)pU* M , t ^)] = Tr[a s [4 )t pC7* t ] (120) 

uniformly in x. 

Proof. According to Lemma EOl Equation (|12U|) is equivalent to 

lim Tr{U* Mit a s U M ,t(x)p} = ^IK t a S U x ,p] (121) 

To prove this let us introduce the notations 

X M {x) = {N+pt)-iu* Mtt {x)a s U M ,t{x)i X(x) = (N+pl)-iui t a s U x , t , (122) 
and show first strong convergence of the defined operators. 
Lemma 3.16. For any r\ £ T-L^ we have 

Jim \\(X M (x)-X(x))ri\\=0 (123) 

uniformly in x. 

Proof. We use the same strategy as for Proposition 13 . 1 01 and decompose rj € Hoo as 

m = ]T 77 (n V n , m = E ^ n) ^ IMI < 6 > ( 124 ) 

n<n e n>n c 

cf. formula ([M]). Then 

\\(X m (x)-X(x))t ] \\ < \\(X M (x)-X(x)) m \\ + \\(X M (x)-X(x))r)2\\ (125) 
< £ h (ll) lll(^M(x) - X(o:))Vn|| + 2Ke (126) 



<(E h ( " ) l)!l(N + dl)^||- + 2X £ , (127) 

Tic 



where the second inequality is the consequence of Lemma l3.121 whereas for the third one assumes 
that M > max n<ne (M e ) in the statement of Lemma [3.111 Since (J2 n<n \f]^ n '\) < n e and (N + 
dT)~2 is bounded, the convergence is proved. □ 

The sequence Xm converges strongly and is norm bounded. This implies for any trace class 
operator convergence of traces Ti(Xm(x)p) with limj\/^oo Tr(X nj(x)p) ~ Tr(X(x)p). Hence we 
get 

Jim \Tr(X M {x)p) - Tr(X(x)p) \ = (128) 
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uniformly in x. Therefore we get with a Schwartz operator p 

Hm Tz{U* Mt {x)a S U M> t{x) p)= lim Tr (( N + P 1 ) 1 ( N + P 1 )~ i U* M t (x) a s U M ,t(x) p) (129) 

M —>oc M— 5-oo 

= lim Tr (X M pCN+pt)i) (130) 
= lim TrX M p = TrXp = Tr(C/* t a s J7 2:it/ 9) (131) 

Af— ^oo 

where p is again a Schwartz operator (and therefore trace class) according to Proposition l3.ll □ 

3.4. Proof of the main theorem 

Based on the result derived so far we can now start to prove the main theorem. This means in 
particular that we have to take the x-dependence stronger into account. We start with a short 
lemma which will be the main tool throughout this Subsection. 

Lemma 3.17. For each M £ N U {oo}, x £ [0,1], t £ M. with \t\ sufficiently small and each 
multiindex S £ { — 1, l} d (including the trivial case \S\ = d = 0) let us define 

u S m,Up) = ^{U* M , t {x) a s U M , t (x) p) (132) 

and consider S C S{T-L) bounded, i.e., 

WpiN+piy/% <K p Vp£X. (133) 

Then we have 

\"M,t, x (p)\ < KK p (134) 

with the constant K from Lemma \3.1'£\ and p £ N sufficiently big ( such that the statement from 
Lemma \3.12\ holds). Moreover ujfj t x (p) is continuous as a function of x. 

Proof. Since p is a Schwartz operator we can rewrite cjf^ t x as 

«&,*,*(/>) = Tr((N+pl)^/ 2 C/ Mt (^)*a 5 C/A/, t (x)p(N + pi)"' 2 ) (135) 

hence 

\<*M,t,*(p)\ < +pl)- p/2 UMAxya S U M Ax)\\\\p(N +pl) p/2 \\i <K-K p (136) 

The first term is bounded by K by Lemma l3.121 the second term is bounded by K p by assumption 
(where p has to be chosen large enough such that Lemma 13.121 is true). This shows boundedness. 

To show continuity in x consider first the case M = oo. For a quadratic Hamiltonian 
cxp(—itH)a s cxp(itH) is a linear differential operator in x with coefficients smoothly depending 
on t. Hence the continuity of x t (p) in x follows from that of exp(—ixtH) for each t fixed. 

If M is finite the situation is effectively finite dimensional, since UM,t( x ) ac ^ s as the identity 
on all ip n with n > Mx; the monomial a s can map ip n to hence, we have to look at the 

Hilbert space spanned by tp n , n = 0, 1, 2, . . . , xM+\S\. Now, continuity of cuf f x t (p) follows from 
the continuity of the map x > Hm(x) in norm for all fixed A, M, since this continuity implies 
norm continuity of UM,t{ x ) m x t. □ 

Lemma 3.18. Consider a sequence of continuous maps: Rm,\ ■ [0, 1] — > S{%), M £ N such that 

1- {R M ,\{x)\x £ [0, 1],M £ N} is a bounded subset of S(H) 
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2- Rm,\(x) converges inS(H) and uniformly on a neighbourhood I C [0, 1] of X to a continuous 
function R^ : [0, 1] — > 5('H o)- 

Then the sequence 

Tr (U* Mjt (x) a s U M>t {x) R M ,x(x) - 17* t a s U^R^x)) (137) 
is uniformly bounded oniG [0, 1] and converges uniformly to on x G /. 

Proof. Convergence of the Rm,\(x) implies that for each p€ N and each e we can find an M P:<i 
such that 

\\(RmAx)~Roo(x))(^+p1Y /2 \\ < VM>M W (138) 
holds. With Lemma T3. 171 this implies 

\wit t x (Rm,x{x) - Roo(x))\ < K^- = e, Vx G J,VM > M t (139) 

Hence, we write 

\Tr(U* M j(x) a s U M ,t{x) R M ,\{x) - U* >t a s U^R^xj) \ < 

\"M,t,x (RmAx) -Roo(x))\ + \wft, t ,x {Roc (X)) - it>x (Roo (x))\ , (140) 

where the first term was just shown to be < e and the second term converges to zero by Propo- 
sition Note, that both statement are independent of x G I. 

Boundedness follows similarly, if use condition [1] and the fact that the map R 0o (x) G S("Hoo) 
is bounded as well (due to continuity and compactness of [0, 1]). □ 

Lemma 3.19. Consider M G N U {oo}, x G [0, 1], t G R with tj sufficiently small, a multiindex 
S G { — 1, l} d and a continuous map R : [0, 1] — » S(Hoo)- Then the functions 

[0,1] 3i^< t ,(N eC (141) 

are continuous. 

Proof. Consider a sequence Xj G [0, 1], j G N converging to x G [0, 1]. Continuity of R implies 
that for each e > and each r G N there is an j r ,e such that 

\\(N + rl) r p\\ < e (142) 

holds for all p in the set 

£ r , e = {R(Xj) - R(x) | j > i e , r }. (143) 

Hence j > j r!e implies 

||WM,x„ t (^))-^M,x,tW)| < 

ll^.tfcO -^))| + l^.t W)!- ( 144 ) 

The first term on the right hand side is bounded by Ke due to (|142p and Lemma 13.171 The 
second term can be made arbitrarily small due to continuity of uifj x t (p) in x; cf. again Lemma 
I3TT71 □ 
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Lemma 3.20. The sequence 

M h> {a M (x) s - x^a s ){N + 2|S|I) -|S| . (145) 

converges to zero in norm uniformly in x € [0, 1] any closed subinterval of x G / C [0, 1] not 
containing 0. 

Proof. The statement is proved if we can show that for each e > we can find an M £ which is 
independent of n and x such that 

\\a s M (x) - x^a s )(N + 2|S|1)-I S ^„|| < ee (146) 

holds for all M > M e . Using the definition of bm(i) in Equations (|38p and (|39p and properties 
of the standard creation and annihilation operators it is easy to see that the elements of the 
sequence in (|146[) are products of l^l terms of the form 



(PmAx, n + r)- y/x)Vn~+P{n + 2 1 ^ | ) - 1 (147) 

which differ only in the parameter r, r'\ we have used a similar argument already in the proof of 
Lemma |3~81 cf. Equation (|84p . The only difference are l^l terms (n + 2|S'|) _1 we have distributed 
among the S factors. Using the bound in lemma [X71 we get 



/3 M ,\(x,n + r) - y/x ^Jn + r' 1 n + r n + r' 

^n + 2\S\ y/n + 2\S\ ~ VM V « + 2\S\ n + 2\S\ 

The right hand side can be bound in an n-independent way to 



n + r n + r' 



< N /max(l, r/2S) max(l, r'/2S). (149) 



y n + 2\S\ n + 2|S| 

This leads to the estimate in (|146[) which concludes the proof. □ 

Lemma 3.21. For each p £ N, t € K sufficiently small and a uniformly bounded sequence 
Rm ■ [0, 1] — > S(Hoc) the sequence 

||(N + 2ptfU M , t {x)R M {x)U M ^x)\\ l Men (150) 

is uniformly bounded in x. 

Proof. Since (N + 2pT) p is a polynomial in a and a* it is sufficient to look at 

||a s C/ M ,t(x)i? A/ (x)[/ A * /jt (x)||i. (151) 

In addition we can rewrite this expression as 

Wa^M^RMixMt^x)^ = \\Ui Itt {x)a s U M Ax)RM(x)\\i (152) 

= ||C/^ t (x)a s [/ M ,t(a ; )(N+p'l)-f'/2 (N+ y I)P 72 i?M(x) || 1 (153) 

< ||C/^ t (x) a s C/ M ,t(a ; )(N+p'l)-f'/2|||| (N + y I)P 72 i?M(:E) || 1 < K 'K X (154) 

The first factor is bounded according to Lemma l3.12l thc second because Rm is uniformly bounded 
by assumption and p *-> a s p is a continuous map on S(Hoo); cf. Proposition 13. II □ 
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Lemma 3.22. For each M 6 N, t 6 K sufficiently small and R G { — 1, l} d the function 

[0, |Tr(C/ M , t (x)i? A/ ( a; )f/; /it (x)(af f ( :E ) - £*«,(;»)*)) | £ C (155) 

is continuous. 

Proof. We can write this expression as a linear combination of therms oj m x t (Rm (x) J with 
coefficients depending continuously on sc. Hence the statement follows from Lemma 13.191 □ 

Now we are ready to finish the proof of Theorem 12.101 To this end let us go back to Equation 
([53)) which was 

\Ti(U M ,tPMU* Mtt a^) - Tr(U x , tPoo Ul t a R )\ < 

) R )) \nM{dx)+ 

Jo 

[ |n((E^, t iMx)Z7? t - ^M,t(a;)i2oo(a;)^ >t (a;))a 00 ( : r) R )|/i M (da:) (156) 
Jo 

where aoo(x) = x' R " 2 a according to Equation (j5Tj) . 

The integrand of the second integral is - according to Lemmas 13.181 and 13.191 - continuous and 
uniformly bounded in x £ [0, 1], and it converges uniformly on x G I to 0. In addition we know 
that the measures [1m converge weakly to the point measure at A. Hence the second integral 
vanishes in the limit M —> 0. 

The integrand of the first integral is continuous as well (Lemma I3.22[) and we can rewrite it 
according to 

\Tr((a M (x) - x R / 2 a R )(N + 2|S|1)-I S I (N + 2\S\1)\ s \U m MRm{x)U m ^x)) \ 

< \\(a R (x) - x R / 2 a R )(N + 2\S\1)~^ \\ \\ (N + 2|S|l)l s l[/ M ,t(a;)i?M(x)f/^ t (x)|| 1 . (157) 

The first factor on the right hand side converges to uniformly in x (Lemma I3.20p while the 
second is uniformly bounded (Lemma 13. 2 1|) . Hence in the limit M — > oo this integral vanishes as 
well, which concludes the proof. 

4. Outlook 

We have shown that a one degree of freedom, continuous variable quantum system can be simu- 
lated by mean field fluctuations of an ensemble of qubits. Our results substantially exceed existing 
schemes. This includes in particular the range of states we can simulate and the treatment of 
the dynamics. Nevertheless, there are a number of open questions which needs to be treated in 
forthcoming papers. 

• Dynamics. Our results concerning the dynamics have two essential restrictions: The state- 
ments are restricted to quadratic Hamiltonians and to small times. The latter is probably 
only a restriction of the methods in the proof and can be removed with a more careful anal- 
ysis of the limiting behavior of the unitaries UM,t(x). One possibility is to show that the 
time evolved density operators Um ,t(x) Rm (x)Ulj t (x) converge in the same as the §Rm(x) 
do, because in this case we can stack arbitrarily many finite time intervals together. 

The discussion of more general Hamiltonians is most likely a more difficult problem. In 
Scction [3.2l we have made heavy use of the fact that Hermite functions are analytic vectors 
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for all quadratic Hamiltonians, and this argument has to be replaced somehow. A mini- 
mal requirement is in any case that the dynamics leaves the space of Schwartz functions 
invariant. The questions arc: Which Hamiltonians have this property and is it sufficient to 
prove an analog of Theorem 12. f Ok ' 

The reference state. The reference state 9 plays a very crucial role in the theory, which 
was overlooked in previous publications. It provides the parameter A which plays in some 
respect the role of an effective h, although this is not universally true, as we have seen 
during the discussion of Theorem 12.91 In particular the limit A — > seems to be boring, 
because in our scheme the observables Qoo and just scale down with A -1 / 2 and therefore 
all expectation values vanishes if A approaches 0. However, this is not necessarily the end 
of the story. Instead of using sequences of states converging according to Definition 12.81 
we can use more divergent ones such that expectation values do not vanish in the limit 
M — > oo. Alternatively we can rescale Qm{x) and Pm(x) by x~ x l 2 and x 1 ^ 2 respectively. 

Another important aspect arises in connection dynamics. The time evolution we have looked 
at in this paper leave the references state invariant (at least asymptotically in the limit 
M — > oo). It would be interesting to break this property and consider situations where the 
reference state is time-dependent. 

Classical observables. Throughout our analysis we have basically ignored the fluctuation 
operator Fm{ct^). A short calculation along the lines of Section |2~T1 shows that Fm{ct^) 
commutes in the limit M —¥ oo with Qm and Pm, hence it describes a classical observables. 
The fact that Tr((73#) = holds indicates that it is related to A. 

Higher dimensions. To extend of our schemes to systems with arbitrary, finite dimen- 
sional one-particle spaces is a very natural task, and the generalization of the discussion in 
Section [2. II is straightforward. The analysis of permutation invariant states along the lines 
of Section ^. 3l is not. This concerns in particular the embedding of irreducible representation 
spaces of unitary groups into Hilbert spaces L 2 (R d ). 

Gaussian states. A lot of theory is available for Gaussian states and it is not completely 
clear how it is related to our approach. In [Sj it is shown that a translation invariant, 
exponentially clustering state u> of an infinite spin chain leads to Gaussian fluctuations if 
restrictions of ui to finite parts of the chain are used as the sequence pm- In this case it clear 
already from the analysis of correlation functions (as in Section 12. ip that the operators 
Qcx; and P^ can be realized as canonical position and momentum in the Schrodingcr 
representation (modulo factors of the form X^ 1 / 2 ). Also the density operator p^ is easily 
calculated since only one and two point correlation functions are needed. It is, however, 
not clear whether these sequences converge as described in Section 12.31 

Another interesting aspect concerns the relation between Gaussianity and correlations. The 
analysis in [5] uses in a crucial way that the pm arises from restrictions of a thermodynamic 
limit state. Our scheme, on the other hand, does not need this assumption and it would 
be interesting up to which degree absence of long range correlations (in a to be specified 
limiting sense) still implies Gaussianity. Even more interesting (and more ambitious as 
well) is question whether we can distinguish classical correlations and entanglement by 
properties of the limit state p^ . 
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A. Fluctuation operators and representations of tensor algebras 

To prove Proposition 12 . 21 let us consider the tensor algebra A over the vector space V = C 2 : 

A = © V ® M ' V ®° = c - ( 158 ) 

MeN 

Together with the tensor product and the involution 

{x {1) <g> • • • ® x (k) )* = x (k) ® ■ ■ ■ ® (159) 

(where x denotes complex conjugation in the canonical basis of V = C 2 ) the space A becomes a 
*-algebra. The vector space V is naturally embedded in A and its canonical basis 

is a complete set of generators. 

The fluctuation operators Fm give now rise to a *-representation of A by 

®M(q)=Q M , $m(p)=Pm (161) 

where Qm = 2~ 1 / 2 -F'm(o'i), Pm = 2~ 1 ^ 2 Fm(o'2) are the operators defined in Equation Any 
density operator pu e S(H® M ) leads to a state 

Wm{X)=Ti{§m{X)pm\ XeA (162) 

of the algebra and if the sequence pm, M £ N has vM fluctuations the limit 

W(X) lim W M (X) = Tr($ M (X)p M ) (163) 

M— yoo 

exists for all X and defines again a state W of A. 

The state W defines the GNS representation^ {%Wi Dw, ^w, &w) of A consisting of a Hilbert 
space Hw, a dense subspace D\y, a *-morphism ttw into the algebra of (unbounded) operators 
Dw — > Dw and a cyclic vector Q,y/ satisfying 

{n w ,7r w (X)Q w } = W{X). (164) 

Without loss of generality we will assume in the following that Dw coincides with the subspace 
generated by vectors ttw(X)Qw £ Dw ■ Applying ttw to the generators q,p from Equation (|lfciU|) 
we get two (unbounded) operators 

Qoo = nw(q), Poo = nwip) (165) 



5 Usually the GNS representation is only considered for C*-algebras, but it can easily be generalized to more 
general *-algebras. We only have to allow representations in terms of unbounded operators. Our case is basically 
an instance of the celebrated reconstruction theorem of Wightman quantum field theory 1 1 21 . 
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with domain D\y . For any polynomial f(q,p) they satisfy 

lim Tr(f(Q M ,P M )pM) =Tr(f(Q 00l P 00 )p 0C ) with p^ = \Cl w )(Sl w \. (166) 

M — too 

This establishes all statements in Proposition 12.21 except the last commutation relations (item 

u. 

To show the latter, let us consider two elements X\, X 2 G A and the expression (cf. Equation 
©) 

Tt(<P m (X 1 )[Qm,Pm]^m(X 2 )) - iXTr(^ M (X 1 )^ M (X 2 )p M ) 

Tr(^ M (X 1 )F u (a 3 )^ M (X 2 )). (167) 



2VM 

If the sequence pu has y/~M fluctuations we have 



lim Ti-(^ M (X 1 )F M (a 3 )'i>M(X 2 )) = C < 00, (168) 



hence 



W(X 1 {qp-pq-i\l)X 2 ) 

= lim [Tr(&M(X 1 )[Q M ,P M ]*M(X2))-i\Tr{*M(X 1 )* M (X 2 )pM)]=0- (169) 

In other words, the limiting state W vanishes on the ideal generated by qp~pq — iXl which leads 
to 

[<2oo,Poo]0 = iA(/> V0GL>w (170) 

as stated. 
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